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ABSTRACT 

The  work  reported  in  this  paper  is  a  part  of  on-going  studies  to  clarify  how  and  to  what  extent  soil  electromagnetic 
properties  affect  the  performance  of  induction  metal  detectors  widely  used  in  humanitarian  demining.  This  paper 
studies  the  specific  case  of  the  time-domain  response  of  a  small  metallic  sphere  buried  in  a  non-conducting  soil 
half-space  with  frequency-dependent  complex  magnetic  susceptibility.  The  sphere  is  chosen  as  a  simple  prototype 
for  the  small  metal  parts  in  low-metal  landmines,  while  soil  with  dispersive  magnetic  susceptibility  is  a  good 
model  for  some  soils  that  are  known  to  adversely  affect  the  performance  of  metal  detectors.  The  included  analysis 
and  computations  extend  previous  work  which  has  been  done  mostly  in  the  frequency  domain.  Approximate 
theoretical  expressions  for  weakly  magnetic  soils  are  found  to  fit  the  experimental  data  very  well,  which  allowed 
the  estimation  of  soil  model  parameters,  albeit  in  an  ad  hoc  manner.  Soil  signal  is  found  to  exceed  target  signal 
(due  to  an  aluminum  sphere  of  radius  0.0127  m)  in  many  cases,  even  for  the  weakly  magnetic  Cambodian  laterite 
used  in  the  experiments.  How  deep  a  buried  target  is  detected  depends  on  many  other  factors  in  addition  to  the 
relative  strength  of  soil  and  target  signals.  A  general  statement  cannot  thus  be  made  regarding  detectability  of 
a  target  in  soil  based  on  the  presented  results.  However,  computational  results  complemented  with  experimental 
data  extend  the  understanding  of  the  effect  that  soil  has  on  metal  detectors. 

Keywords:  landmine  detection,  humanitarian  demining,  electromagnetic  properties  of  soil,  magnetic  suscepti¬ 
bility  of  soils,  magnetic  viscosity,  pulsed  induction  metal  detectors 


1.  INTRODUCTION 

The  electromagnetic  induction  metal  detector  is  one  of  the  most  commonly  used  tools  to  detect  buried  landmines 
in  humanitarian  demining.  One  of  the  important  factors  that  can  influence  the  signal  produced  in  such  a  detector 
is  the  host  soil.  A  systematic  analytical  framework,  based  on  well-established  techniques  in  geophysics  and  non¬ 
destructive  testing,  for  studying  the  effects  of  soil  electromagnetic  properties  was  presented  in  an  earlier  paper  [1], 
including  a  discussion  of  the  relevance  and  requirements  of  soil-related  research  in  humanitarian  demining! 
Results  of  further  study  and  additional  numerical  results  were  included  in  a  subsequent  paper  [2].  Most  of 
the  discussion  in  these  previous  papers  has  been  based  on  frequency  domain  analysis.  Since  a  large  number 
of  detectors  used  in  practice  employ  pulsed  transmitter  currents,  studies  in  the  time  domain  are  considered 
essential.  To  that  end,  the  present  paper,  which  is  a  part  of  on-going  studies  to  clarify  how  and  to  what 
extent  soil  electromagnetic  properties  affect  the  performance  of  metal  detectors,  analyzes  the  specific  case  of  the 
time-domain  response  of  a  small  metallic  sphere  buried  in  a  soil  with  frequency-dependent  complex  magnetic 
susceptibility.  The  sphere  is  chosen  as  a  simple  prototype  for  the  small  metal  parts  in  low-metal  landmines,  while 
soil  with  dispersive  magnetic  susceptibility  is  a  good  model  for  some  soils  [3]  that  are  known  to  adversely  affect 
the  performance  of  metal  detectors. 

The  analytical  results  and  framework  for  computation  are  summarized  in  Section  2,  which  also  includes 
a  discussion  of  the  specific  model  for  soil  magnetic  susceptibility  used  in  this  paper.  A  brief  description  of 
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experimental  set  up  used  to  measure  the  induction  response  of  soil  and  targets  is  given  in  Section  3.  Numerical  and 
experimental  results  for  soil  and  target  response  are  discussed  in  Section  4.  A  summary  and  overall  conclusions 
are  given  in  Section  5. 


2.  ANALYSIS 

The  geometry  employed  to  analyze  the  metal  detector  problem  is  the  same  as  the  one  used  in  [1,2]  and  is 
repeated  in  Figure  1.  The  soil  is  modelled  as  a  homogeneous  halfspace  and  the  metal  detector  head  is  taken  to 
consist  of  a  pair  of  concentric  and  coplanar  circular  coils.  The  transmitter  coil  radius  and  number  of  turns  are 
respectively  a  and  Nt  while  those  of  the  receiver  coil  are  b  and  Nr.  The  coordinate  systems  used,  the  magnetic 
permeability  (fij),  electrical  conductivity  {p j )  and  permittivity  tj  of  the  various  regions  are  shown  on  the  figure, 
where  the  subscript  j  =  0, 1, 2  indicates  air,  soil  and  target  material  respectively.  The  coil  assembly  is  at  a  height 
z  —  h  from  the  air-soil  interface.  In  the  frequency  domain,  the  transmitter  current  is  defined  to  be  It with 
u)  =  27r/,  where  /  is  the  frequency  of  operation.  Time  domain  results  will  be  computed  from  known  frequency- 
domain  ones  through  inverse  transforms.  Except  for  some  simple  cases,  where  analytical  results  are  available, 
inversions  will  be  carried  out  using  a  numerical  technique  presented  in  [4]. 


Object  (^ ,  a2  ,  ^ ) 


Figure  1.  Geometry  of  the  problem. 


2.1.  Soil  Model 

In  this  paper,  only  the  specific  case  of  non-conducting  soil  (<xi  =  0  S/m)  with  frequency-dependent  magnetic  sus¬ 
ceptibility  *(cj)  will  be  considered.  As  well,  the  permittivity  of  soil  will  be  assumed  to  be  frequency-independent, 
the  same  as  that  of  free  space  and  will  be  not  be  considered  in  the  results  (i.e.,  displacement  currents  will  be 
ignored).  The  particular  model  of  magnetic  behaviour  used  in  the  analysis  is  as  follows: 


Hi  =Mo(l+x(u>)) 

=  x'(u>)  +  =  xo  - 


•  In 


iu)r2  + 


ln(r2/ri)  iu>n  + 


i) 


(1) 
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where  /to  and  m  are  the  magnetic  permeability  of  free-space  and  soil  respectively.  x(w)  is  the  complex  frequency- 
dependent.  susceptibility  of  soil  with  real  and  imaginary  components  x'  and  x"  respectively,  while  xo  is  the  d.c. 
value  of  susceptibdity.  tx  and  r2  are  the  lower  and  upper  limits  respectively  of  magnetic  relaxation  time  constants 
of  the  soil.  This  model,  which  has  been  discussed  in  previous  work  [1, 2, 5-9],  is  often  used  to  represent  soils  that 
display  magnetic  viscosity.  Such  soils  are  common  in  many  landmine-affected  regions  of  the  world. 

Measurements  allowing  unique  determination  of  all  three  parameters  of  the  model,  namely,  xo,  Ti  and  r2 
are  not  readily  available.  (Nevertheless,  an  ad  hoc  trial-and-error  technique  of  estimating  xo,  Ti  and  r2  from 
available  frequency  and  time-domain  data  is  discussed  in  Section  4.1).  However,  an  approximation  derived  from 
Equation  2  can  often  be  used  in  conjunction  with  measurements  made  by  commercial-off-the-shelf  equipment  to 
obtain  results  that  are  useful  in  practice.  If  measurement  times  and  parameter  values  are  such  that  wti  «  1 
and  wr2  »  1  (or  equivalently  t/ n  »  1  and  t/r2  «  1),  it  can  be  shown  from  Equation  2  that 

W  _2  xo 

d  Ituj  7r  /n(T2/ri)  (3) 

According  to  this  relationship  which  has  been  discussed  in  detail  by  others  [8, 10, 11],  the  slope  (with  respect 
to  lnu>)  of  the  real  part  of  susceptibility  is  constant  and  is  proportional  to  its  imaginary  part.  In  practice, 
measurements  of  x'  at  two  frequencies  provided  by  an  instrument  such  as  the  popular  Bartington  MS2B  dual¬ 
frequency  susceptibility  meter  [12],  which  measures  x'  at  465Hz  and  at  4650Hz,  are  used  to  estimate  the  required 
slope  assuming  linear  variation  with  respect  to  lnu>.  This  approximation  will  be  used  in  parts  of  this  paper 
to  calculate  time-domain  response  of  soil.  It  should  be  noted  that  measurements  of  complex  susceptibility 
over  a  much  wider  frequency  band  [11, 13],  albeit  for  a  limited  number  of  soil  samples,  indicate  that  the  linear 
approximation  may  be  adequate  for  some  real-world  soils. 


2.2.  Response  of  Soil 

Using  results  from  [1,2]  one  can  write  down  the  expressions  for  induced  voltages  under  various  conditions.  The 
voltage,  normalized  by  the  product  NtNrI,  induced  in  the  receiver  coil  by  a  non-conducting  (cr,  =  0)  but 
magnetic  halfspace  is  given  by 


vsml(Lj)  =  i/i0u>nab 


(4) 


where  m{h)  —  fQ  Ji(Aa)Jr1(A6)exp(-2A/i)dA,  J\  denotes  a  Bessel  function  of  the  first  kind  and  order  1,  and  A 
is  an  integration  variable.  Using  integral  no.17  on  page  316  in  section  13.4.1  of  [14]  m(h)  can  be  shown  to  be 


(5) 


where  k2  =  (a+$b+4hi ,  and  K  and  E  are  complete  elliptic  integrals  of  the  first  and  second  kind  respectively, 
defined  as  K  =  J”/2  (l  -  k2sin20)~1/2 d0  and  E  =  f*/2  (l  -  k2sin2$) 1/2 d0. 

The  frequency-domain  expression,  Equation  4,  can  be  used  to  derive  time-domain  relations  through  the  use 
of  Laplace  and  Fourier  transform  theory.  For  example,  the  response  to  a  step-current  excitation  can  be  derived 
by  calculating  the  inverse  Fourier  transform  of  (u) / (iw)  or,  equivalently,  by  calculating  the  inverse  Laplace 

transform  of  vsm  (s)/s  where  the  transform  variable  is  defined*  by  s  =  iw.  A  numerical  implementation  of 
this  approach  using  limited  frequency-domain  data  is  given  in  [4].  The  author  is  not  aware  of  any  analytical 
expression  for  the  required  inverse  transforms  when  no  approximation  is  made  regarding  x(w)  in  Equation  4  and 
m  such  cases  use  will  be  made  of  the  numerical  technique  [4]  already  noted.  However,  if  the  soil  is  very  weakly 

magnetic  (  |xl  «  1)  which  is  commonly  the  case  with  natural  soils,  the  response  for  t  >  0  to  a  step-current 
excitation  that  at  t  =  0  can  be  written  as  [1] 


«££  (t)  =  Honabm(h)—^—  .  ± 

_ _  21n(r2/r1)  t 

'Strictly  speaking  v,,oi,(u)  and  x(w)  should  be  written  as  vsoil(iu)  and  x(M  respective 


(6) 
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when  measurement  times  and  parameter  values  are  such  that  ^  >>  1  and  ^-  <<  1.  The  discussion  in  the 
literature  often  centers  around  the  step  response.  However,  the  response  due  to  finite  current  pulses  which  are 
used  in  practice  must  be  considered  in  the  analysis  of  metal  detectors.  One  common  current  waveform,  which 
is  also  the  shape  used  in  the  experiments  for  this  work,  is  a  short  triangular  pulse  of  duration  T.  Starting  from 
the  above  step  response  (Equation  6)  and  applying  Laplace  transform  theory!,  one  can  show  that  the  response 
to  a  triangular  current  pulse  of  duration  T  and  peak  value  70  is  given  by 


soil 

ramp 


«)  =  .  h  (,„! ±1  - 


21n(r2/ri) 

where  t  is  measured  from  current  turn-off. 

Further  simplifications  can  be  obtained  by  using  Equation  3  in  Equations  6  and  7: 

Ktcp  (t)  =  -Ho*abm(h)—  ■  7 
7 r  t 


(7) 


=  —^fiQTrabm(h)  • 

dx' 

d  Inuj 

1 
*  t 

vr7mp  (<)  =  -Honabm(h)  — 

7T 

’■■H 

In— 

=  -““Mott abm(h)  ■ 

dx' 

d  Inu) 

lo 
'  T 

(8) 


-?) 


(9) 


2.3.  Response  of  Sphere 

In  this  paper  only  the  case  of  a  non-magnetic  conducting  sphere  lying  on  the  axis  of  the  coils  will  be  considered. 
Following  [1],  the  voltage,  normalized  by  the  product  NtNrI,  induced  in  the  receiver  coil  by  such  a  target  buried 
in  non-conducting  soil  (<T\  =0,  which  makes  the  equivalent  current  dipole  moment  P4,  =  0)  can  be  written  as: 

«sphcrcM  =  iuJHoHblzMz  (10) 

where  H\z  is  the  *  component  of  magnetic  field  intensity  that  would  be  produced  in  soil  at  the  location  of 
the  sphere  centre  if  a  unit  current  were  flowing  in  the  receiver  coil  of  radius,  b.  Mz  is  the  z  component  of 
the  scatterer  equivalent  magnetic  dipole  moment,  which  for  a  sphere  of  radius,  R,  and  electrical  parameters, 
(M2>C2,<r2)  embedded  in  a  half-space  with  parameters  (/*i,  £1,01)  can  be  written  as: 


M  jj..  x  ^  7^  )  ~i~  2/x2]  sinh^lfyR)  (2/^2  -f  cosh^tc^R) 

[Mi(l  +  fc22 R2)  —  fJ.2]  sinhfaR)  +  (^2  —  /ix)/c2I?  cosh(fc2J?) 

where  Hlz  is  the  field  produced  at  the  sphere  centre  by  unit  transmitter  current,  and  k22  =  ia2p2u>.  The 
expression  for  the  in-soil  field  H\z  is  given  in  [1,2],  The  field  H\z  can  be  also  obtained  from  the  same  expression 
by  simply  replacing  a  with  b  and  setting  7=1. 

As  in  the  case  of  soil,  time-domain  responses  can  be  calculated  by  applying  inverse  Fourier/Laplace  transforms 
to  appropriate  frequency  domain  expressions.  For  example,  inverse  Laplace  transform  oivslpheTe(s)/s  would  result 
in  the  response  of  the  sphere  to  a  step-current  excitation.  The  author  is  not  aware  of  an  analytical  solution  when 
the  sphere  is  buried  in  a  soil  with  general  magnetic  properties  characterized  by  Equations  1  and  2.  In  such  cases, 
inverse  transforms  will  be  calculated  numerically.  However,  if  the  sphere  is  situated  in  free  space,  analytical 
solutions  are  possible.  A  detailed  theoretical  and  experimental  study  of  the  time  domain  response  of  a  sphere, 
which  considers  both  permeable  and  nonpermeable  spheres,  higher  order  induced  multipoles  (order  of  spatial 
mode  n  >  1),  the  effect  of  pulse  shape  and  measurement  electronics  with  coils  of  finite  thickness,  is  given  in  [15], 

fAn  alternate  derivation  is  discussed  in  [9]. 
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Here  only  the  simpler  case  which  considers  a  nonpermeable  sphere  (fi2  =  (M>)  and  only  induced  dipoles  (n  =  1) 
and  assumes  coils  to  have  no  thickness  will  be  treated.  Following  the  technique  discussed  in  [15],  the  step  response 
or  the  sphere  in  free  space  can  be  written  as 


sphere  /,\  _  3nR3(ioa2b2  1  n,y  f 

Vstcp  ~ W+^WW+W^ "p2^6  9 

m=  1 


.where  02  =  ii0a2R2 


(12) 


Similarly,  one  can  show  that  the  response  to  a  triangular  current  pulse  of  duration  T  that  drops  to  0  from 
its  peak  value  of  I0  at  t  =  0  is  given  by 


.  sphere  (f\  _  fr> 
uramp  \L )  ~  TpT 


3nR3Hoa2b2 


T  ( h 2  +  a2)3/2(/i2  +  62)3/2 


-£  +  n+±  -  £  t  .-^«) 


(13) 


where  t  >  0  and  U(t)  is  a  unit  step  function.  Although  there  would  be  no  occasion  to  use  this  particular  result 
m  this  paper,  it  is  included  here  for  completeness  and  to  indicate  that  one  should  be  mindful  of  the  pulse  shape 
when  analyzing  metal  detectors.  Future  work  will  use  this  result. 

3.  DESCRIPTION  OF  EXPERIMENTS 

A  report  incorporating  details  of  the  experimental  setup  is  being  prepared.  Only  a  minimal  description,  sufficient 
to  understand  the  parameters  of  the  experimental  data,  is  given  here.  All  measurements  were  conducted  using 
a  suitably  modified  Scheibel  AN19/2  metal  detector.  The  sensor  coils  had  the  following  nominal  parameters 

(  A  nnoSmit*ter  '  20  tUmS  With  a  mean  radius  of  0  1275m;  (2)  receiver  coil  -  33  turns  with  a  mean  radius 
of  0.093m  An  external  function  generator  was  used  to  drive  the  transmitter  to  produce  a  unipolar  triangular 
current  pulse  tram  with  the  following  parameters:  pulse  duration,  T=132  fis;  maximum  current,  70  =3  06  A 
and  a  pulse  period  of  ~150  ms.  The  receiver  amplifier  had  a  d.c.  gain  of  270  and  a  3-dB  bandwidth  of  200  kHz! 
The  receiver  output  was  sampled  every  3  (js  and  digitized  by  a  16-bit  A/D  converter.  To  obtain  improved  S/N 
ratios,  5000  cycles  of  the  output  were  coherently  averaged  to  arrive  at  a  single  response  measurement  The 
transmitter  current  waveform  was  simultaneously  measured  in  a  similar  way.  The  background  response,  with  the 
sensor  head  positioned  away  from  the  ground  (at  a  height  at  least  1.5  m  above  the  ground)  and  other  objects 
was  periodically  measured.  The  background  response  was  subtracted  from  the  total  response  in  the  presence  of 
an  object  (e.g.,  soil,  metal  target)  in  order  to  obtain  the  desired  secondary  response  due  only  to  the  object(s). 
for  the  data  reported  in  this  paper,  the  target  used  was  an  aluminum  sphere  of  radius  R  =  0  0127  m  and 
conductivity  a2  =  3.54  x  107  S/m. 

In-air  target  response  was  measured  in  a  climate-controlled  laboratory  of  nonmetallic  construction.  Soil 
response  and  m-soil  target  response  were  measured  using  a  2  m  deep  soil  pit  with  a  4  mx6  m  surface,  filled 
with  Cambodian  laterite.  Susceptibility  of  soil  samples  was  measured  using  the  Bartington  MS2B  sensor  at  two 
frequencies  465  Hz  and  4650  Hz.  Six  soil  samples  from  the  patch  of  soil  over  which  measurements  were  made 
were  collected  and  their  susceptibility  measured.  One  of  these  samples  were  taken  from  the  spot  directly  under 
the  center  of  the  coils  and  the  remaining  five  were  all  within  a  transmitter  coil  radius  of  this  point.  The  average 
of  the  susceptibilities  of  these  samples  is  assumed  to  represent  an  “equivalent”  susceptibility  for  the  halfspace 
model.  The  accuracy  of  this  assumption  will  depend  on  the  actual  degree  of  uniformity  of  the  soil  and  on  sensor 
height.  The  measured  susceptibility  values  are  shown  in  Table  1. 


4.  RESULTS 

4.1,  Soil  Response  and  Estimation  of  JVTodel  Parameters 

Measured  and  theoretical  responses  of  the  Cambodian  soil  are  compared  in  Figure  2  for  a  number  of  sensor 
heights  (6).  Numerical  techniques  of  transforming  frequency-domain  results  to  the  time  domain,  when  triangular 
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Table  1.  Measured  real  part  of  magnetic  susceptibility  of  soil  samples  in  SI  units,  a  multiplying  factor  of  1CT5  being 
implicit. _ 


Sample  No. 

Real  x  at  465  Hz,  xlf 

Real  x  at  4650  Hz,  xhf 

Xlf  -  Xhf 

1 

252 

238 

14 

2 

208 

192 

16 

3 

320 

306 

14 

4 

249 

'  234 

15 

5 

227 

217 

10 

6 

295 

279 

16 

Average 

258.5 

244.3 

14.2 

waveforms  like  the  ones  used  in  the  experiments  are  used,  are  still  under  development.  Partly  for  that  reason, 
the  theoretical  responses  presented  are  only  the  ones  calculated  using  the  approximate  relation,  Equation  9.  The 
s*°Pe  Tuu3  *s  estimated  from  the  average  values  of  xlf  and  Xhf  shown  in  Table  1,  and  the  overall  gain  of  the 
measurement  system  is  accounted  for  in  calculating  the  response.  Given  the  approximations  involved  in  deriving 
Equation  9,  experimental  uncertainties  and  inhomogeneity  of  soil  properties,  the  agreement  between  measured 
and  calculated  responses  is  considered  remarkable.  The  data  also  show  the  rapid  decrease  of  soil  signal  with 
sensor  height,  particularly  for  times  close  to  current  turn-off. 

The  good  agreement  seen  above  means  that,  when  the  sought-after  response  can  be  cast  as  a  function  of 
or  x"  as  in  Equation  9,  one  need  not  explicitly  know  the  values  of  the  individual  soil  model  parameters,  namely, 
Xo,  T\  and  r2.  However,  in  general,  knowledge  of  the  values  of  the  individual  soil  parameters  may  be  needed. 
For  example,  the  equations  used  in  this  paper  for  calculating  fields  inside  the  soil  halfepace  and  the  response  of  a 
buried  target  require  knowledge  of  the  individual  parameters.  These  equations  will  be  modified  in  the  future  to 
take  advantage  of  the  approximate  linear  (with  respect  to  Irud)  behaviour  of  \  m  some  natural  soils,  as  is  done 
in  deriving  Equation  9.  For  now,  the  general  formulation  will  be  used. 

The  author  is  not  aware  of  any  reliable  and  accurate  method  of  estimating  the  individual  soil  parameters 
from  measurements  provided  by  a  readily  available  instrument.  Instruments  being  developed  [13]  to  measure 
complex  susceptibility  of  soil  over  a  wide  frequency  spectrum  should  allow  estimation  of  individual  parameters 
in  the  future.  In  this  paper,  an  ad  hoc  method  is  used  to  get  a  workable  estimate  of  the  parameters  through  a 
trial- arid-error  process.  The  method  followed  is  described  below. 

A  measured  response  (the  measurement  at  h  =  Om  was  used)  is  linearly  fit  to  the  function  represented  by 
Equation  7  to  estimate  the  factor  .  The  fit  obtained  is  very  good  and  the  uncertainty  in  the  estimated 

value  of  ln(*/Tl)  small.  Obviously,  knowing  an  accurate  value  for  this  factor  does  not  allow  unique  determination 
of  the  individual  parameters.  Knowing  the  general  shape  of  frequency  variation  of  susceptibility,  a  value  of  xo 
is  arbitrarily  chosen  to  be  one  slightly  higher  than  xlf,  in  order  to  restrict  the  range  of  possible  values  of  the 
parameters.  This  choice  fixes  the  value  of  From  this,  various  pairs  of  values  of  r2  and  T\  can  be  selected 
for  further  analysis.  The  values  of  Xo,  Ti  and  r2  thus  selected  are  then  used  as  initial  guesses  in  a  non-linear 
fitting  program  that  fits  Equation  7,  treating  all  three  parameters  as  unknowns,  to  the  same  measured  data  in 
an  attempt  to  estimate  parameter  values.  With  some  trial  and  error  a  good  fit  is  obtained  but,  as  expected,  the 
uncertainty  in  parameters  are  extremely  large.  The  recovered  parameters  are  used  in  Equation  2  to  compute 
X^w),  and  the  calculated  values  are  compared  with  the  two  available  measured  values,  namely,  xlf  and  Xhf- 
The  parameters  are  then  varied  in  a  trial-and-error  manner  to  get  a  good  agreement  with  the  measured  slope 
and  values.  One  set  of  parameters  that  produced  a  reasonable  agreement  for  the  average  values  of  xlf  and  xhf 
in  the  case  of  the  Cambodian  soil  are  Xo  =  0.003,  7i  =  10~22  s  and  t2  =  0.5  s.  This  set  of  parameters  are  simply 
used,  as  needed,  for  the  computations  in  this  paper  and  no  claim  is  made  as  to  their  uniqueness  or  physical 
significance. 

In  the  future,  a  better  method  will  be  developed  to  replace  the  above  somewhat  ad  hoc,  crude  and  cumber¬ 
some  process.  One  of  the  reasons  that  the  non-linear  fitting  process  does  not  result  in  more  certain  parameter 
estimation  is  that  Equation  7  incorporates  the  assumption  of  linearity  (with  respect  lncj)  of  susceptibility  over 
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the  entire  frequency  spectrum.  Equation  7  needs  to  be  replaced  by  an  expression  that  does  not  depend  on  such 
an  approximation  of  Equation  2.  Fitting  such  a  model  to  measurements  may  produce  better  estimates  of  the 
individual  parameters.  Alternatively,  if  most  soils  are  found  from  measurements  to  be  linear  in  behaviour  over 
the  frequency  band  of  interest,  equations  to  compute  various  fields  and  responses  can  be  modified  to  obviate  the 
need  for  individual  values  of  the  parameters.  All  these  aspects  will  be  investigated  in  the  future. 


Figure  2.  Measured  and  calculated  response  of  Cambodian  soil  for  a  number  of  sensor  heights  (h) 
calculated  responses  while  circles  are  measured  data. 


Solid  lines  are 


4.2.  Conducting  Sphere 

As  mentioned  already,  numerical  techniques  of  inversion  of  frequency-domain  results  to  the  time  domain  when 
triangular  waveforms  are  used  are  still  under  development.  For  this  reason,  only  response  due  to  a  step-function 
excitation,  which  can  be  computed  using  the  method  in  [4],  will  be  used  to  study  the  case  of  a  buried  sphere. 
Although  practical  metal  detector  systems,  including  the  one  used  in  the  experiments  here,  do  not  measure  such 
a  response  some  basic  questions  regarding  the  effect  of  soil  can  be  investigated  through  the  use  of  the  ideal  step 
response.  To  this  end,  the  response  of  the  aluminum  sphere  in  air  and  that  when  it  is  buried  in  soil  are  computed 
and  compared  with  the  response  of  the  soil  halfspace  alone.  As  a  check  on  computations,  the  step  response  of  the 
sphere  m  air  is  computed  using  two  methods:  (1)  numerical  inversion  of  vsPherc (w) / (iu)  from  Equation  10;  and 
(2)  direct  computation  using  the  analytical  expression  given  by  Equation  12.  Response  of  the  soil  is  evaluated 
using  numerical  inversion  of  v^(u) /(iw)  from  Equation  4.  Figure  3  shows  the  step  response  of  the  aluminum 
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Figure  3.  Step  response  of  aluminum  sphere  in  air  and  in  Cambodian  lateritic  soil  for  a  number  of  target  depths  d, 
compared  with  the  response  of  the  soil  for  sensor  height  h  =  0  m.  Response  of  the  sphere  in  air  is  computed  by  two 
methods:  (1)  numerical  Fourier  inversion  and  (2)  evaluation  of  analytical  expression.  Soil  parameters  used  are  v0  =  0  003 
n  =  1CT22  s  and  t2  =  0.5  s. 


sphere  in  air  and  in  Cambodian  lateritic  soil  for  a  number  of  target  depths  d  along  with  the  response  of  the 
soil  when  the  sensor  is  lying  on  it,  that  is,  h  =  0  m.  For  the  cases  under  consideration,  an  examination  of  the 
results  shows  the  following: 

(a)  Calculation  of  response  of  the  sphere  in  air  by  the  two  methods  produce  the  same  result,  which  verifies  the 
accuracy  of  the  particular  technique  of  numerical  Fourier  inversion. 

(b)  For  the  soil  in  question,  the  response  of  the  sphere  in  air  and  in  soil  are  not  noticeably  different.  This  would 
allow  the  use  of  the  in-air  response  in  the  place  of  in-soil  response  which  could  be  more  difficult  to  compute 
in  general. 

(c)  The  signal  due  to  the  metal  sphere  dominates  the  soil  signal  when  the  target  depth  is  “shallow”,  d  <  0.10  m. 
But,  even  for  the  soil  with  relatively  low  susceptibility,  the  soil  signal  begins  to  exceed  the  signal  due  to  the 
sphere,  which  has  significant  metal  content  compared  with  that  in  some  low-metal  antipersonnel  landmines, 
when  it  is  buried  deeper  than  about  0.10  m.  At  a  depth  as  low  as  0.10  m,  the  initial  part  of  the  soil  signal 
marginally  exceeds  the  target  signal.  At  d  =  0.15  m  the  soil  signal  completely  dominates  the  target  signal. 
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f  deteCtfab'hty  0f  the  tar8et  can  be  made  just  from  the  above  observations. 

The  output  of  the  detector  is  the  sum  of  the  signal  due  to  the  soil  and  that  due  to  any  target.  If  the  soil  is 

Um  0n7ne  C°U  d’  “  pdnciple’  ^  balance  out  the  signal  due  to  the  soil  and  detect  the  target 
1^17  depth  “  in  air  H™ever-  in  Practice  the  success  of  such  a  simple  approach  will  depend  on  the 
degree  of  spatial  uniformity  of  the  soil,  the  degree  to  which  sensor  standoff  can  be  maintaineda  constant 
the  dynamic  range  of  the  receiver  system,  particular  balancing  and  detection  technique  employed,  and  so 

*1  '  .  'd  fr°m  ,he  qU®St'0n  of  detection,  the  relatively  strong  soil  signal  could  produce  false  detections  if 

there  is  small  scale  spatial  variation  in  soil  susceptibility. 

(e)  The  shape  of  decay  of  soil  signal  is  different  from  that  of  the  target  signal,  a  fact  that  may  be  exploited  for 
discriminating  between  the  two  signals  [3],  at  may  oe  exploited  tor 
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F,gUre  s‘ep  response  of  aluminum  sphere  in  air  and  in  Cambodian  lateritic  soil  for  a  number  of  target  depths  d 

zsr. 1  jssvss;  t  *■ for  *  - 0  03  ■  ■“  *  •  — « — — . 

tvnvTrC  4  PreSe”tS  data/°r  the  same  situations  as  in  Figure  3  expect  that  the  sensor  is  now  used  at  a  more 
typical  operating  height  of  h  =  0.03  m.  Comparison  of  the  two  figures  show  that  most  of  the  comments  made 

sensor  Sdoff  from  the  ^  In  in  FigUre  4’  the  soil  si^al  is  reduced  due  to  the  increased 

h  “d  ff  h  ^  d  SUrfaCC’  Whlch  however  does  not  result  in  a  higher  ratio  of  soil  to  target  signal 

because  there  is  also  an  accompanying  reduction  in  the  target  signal  due  to  the  increased  distance  betweenThe 
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target  and  the  sensor.  In  fact,  for  the  d  —  0.10  m  case,  the  initial  portion  of  the  soil  signal  dominates  the  target 
signal  more  when  h  =  0.03  m  than  when  h  =  0  m.  The  same  comments  about  detectability  as  for  Figure  3  will 
apply  here  as  well. 

The  exact  nature  of  the  effect  of  soil  on  target  signal  will  depend  on  the  soil  parameters  which  have,  so  far, 
been  taken  as  the  ones  estimated  for  the  particular  Cambodian  soil  used  in  the  experiments.  Obviously,  it  is 
not  possible  to  present  results  for  all  possible  values  of  the  parameters  xo,  Ti  and  T2.  To  give  an  idea  of  how 
some  other  values  of  these  parameters  may  affect  the  soil  and  target  signals,  the  case  of  the  sphere  buried  at 
d  =  0.10  m  in  soils  with  four  different  sets  of  parameters  is  shown  in  Figure  5.  The  parameters  used  include 
some  values  used  in  the  author’s  previous  work  [1,2]  when  no  measurements  of  real  soils  were  available.  These 
values  are  used  here  for  the  purposes  of  continuity  and  illustration,  and  no  claim  is  made  as  to  their  occurrence 
in  any  real-world  soil.  Figure  5  shows  the  following: 
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Figure  5.  Step  response  of  aluminum  sphere  in  air  and  in  soils  with  sets  of  parameters  as  follows:  (a)  *o  =  0.05  n  = 

10_6  a,T2  =  °:53  51  (b)  Xo  =  °-5)Tl  =  10"22  9’T*  =  0  5  »;  (<0  Xo  =  0.05,  n  =  10-6  s,t2  =  10-3  s;  (d)  Xo  =  0.5,’n  = 
10  s,T2  —  10  s.  Sphere  depth,  d,  in  all  cases  is  0.10  m,  and  sensor  height,  h,  is  0  m. 

(i)  Because  of  the  much  higher  values  of  Xo  than  in  the  cases  considered  so  far,  the  soil  signal  is  much  higher 
in  general  and  it  exceeds  the  target  signal  by  a  large  amount  for  the  depth  of  burial  considered,  namely, 
d,  =  0.10  m.  The  higher  soil  signal  is  expected  from  Equation  6. 

(n)  For  the  same  value  of  \o,  the  soil  signal  is  higher  for  the  more  lossy  soil,  that  is,  the  one  with  the  lower 
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value  of  t2/ti  (Compare  Figure  5  (a)  with  Figure  5  (c),  and  Figure  5  (b)  with  Figure  5  (d)).  This  behaviour 
is  also  expected  from  Equation  6. 

(iii)  The  effect  of  the  soil  on  the  signal  of  the  sphere  is  very  small  (see  Figure  5(a)  and  (c))  even  for  a  value  of 
Xo  =  0.05,  which  is  more  than  an  order  of  magnitude  higher  than  that  of  the  Cambodian  soil  used  in  the 
experiments.  So,  for  a  lot  of  natural  soils,  one  may  be  able  to  get  away  with  equating  the  target  response  in 
air  with  that  in  soil.  Although  a  xo  of  0.05  is  considered  high  for  natural  soils,  it  is  not  uncommon.  Thus 
if  a  detector  is  able  to  perfectly  cancel  out  the  soil  signal  (not  possible  in  practice),  the  detection  depth  for 
a  target  will  remain  the  same  as  in  air. 

(iv)  For  the  highest  value  of  *o  considered,  namely,  0.5,  the  soil  has  significant  effect  on  the  target  signature 
(see  Figure  5(b)  and  (d)).  The  values  of  relaxation  time  constants  involved  determine  how  the  time-domain 
response  of  the  sphere  is  affected.  It  is  clear  that  for  such  a  soil,  even  if  the  soil  signal  can  be  balanced  out 
perfectly,  the  detection  distance  in  soil  would  be  lower  than  that  in  air. 

A  value  of  Xo  =  0.5  is  of  the  order  of  the  highest  values  found  in  some  “soil  lanes”  used  in  landmine 
detection  research  (for  example,  NVESD  at  Ft.  Belvoir  have  used  a  lane  with  xo  ~  0.35). 

5.  SUMMARY  AND  DISCUSSION 

The  electromagnetic  induction  response  of  a  soil  with  frequency-dependent  magnetic  susceptibility  and  of  a  small 
aluminum  sphere  buried  in  it  were  studied  in  order  to  understand  the  adverse  effect  that  some  soils  occurring  in 
landmine-affected  regions  can  have  on  metal  detectors.  Experimental  data  confirmed  numerical  prediction  that 
real-world  lateritic  soil  (in  this  case  soil  from  Cambodia)  can  produce  a  strong  response  which  can  exceed  that 
due  to  the  buried  metal  target.  Approximate  theoretical  expressions  for  very  weakly  magnetic  soils  were  found 
to  fit  the  experimental  data  very  well,  which  allowed  the  estimation  of  soil  model  parameters,  albeit  in  an  ad 
hoc  manner. 

Although  soil  signal  was  found  to  exceed  target  signal  in  many  cases,  no  definitive  statement  about  the 
detectability  of  a  target  in  soil  can  be  made  just  from  that  observation.  The  output  of  the  detector  is  the  sum 
of  the  signal  due  to  the  soil  and  that  due  to  the  target.  If  the  soil  is  spatially  uniform  one  could,  in  principle, 
simply  balance  out  the  signal  due  to  the  soil  and  detect  the  target  at  the  same  depth  as  in  air,  as  long  as  soil 
susceptibility  is  not  too  high,  which  is  true  for  a  large  number  of  naturally  occurring  soils.  However,  in  practice 
the  success  of  such  a  simple  approach  will  depend  on  the  degree  of  spatial  uniformity  of  the  soil,  the  degree 
to  which  sensor  standoff  can  be  maintained  a  constant,  the  dynamic  range  of  the  receiver  system,  particular 
balancing  and  detection  technique  employed,  and  so  on.  Aside  from  the  question  of  detection,  the  relatively 
strong  soil  signal  could  produce  false  detections  if  there  is  small  scale  spatial  variation  in  soil  susceptibility. 

Future  work  will  include  development  of  a  better  method  of  estimating  soil  model  parameters,  investigation 
of  soil  and  target  response  using  triangular  current  pulse  trains  instead  of  the  step  current,  application  of  the 
analysis  and  experimental  verification  to  more  targets  and  soils,  and  comparison  of  results  with  those  obtained 
by  others  using  different  analysis  techniques. 
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